Quantum phase transitions of 2-d dimerized spin- 1/2 Heisenberg models with spatial 
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Motivated by the unexpected Monte Carlo results as well as the theoretical proposal of a large 
correction to scaling for the critical theory of the 2-d staggered-dimer spin-1/2 Heisenberg model 
on the square lattice, we study the phase transitions induced by dimerization of several dimerized 
quantum Heisenberg models with spatial anisotropy using first principles Monte Carlo method. We 
focus on investigating the finite-size scaling of the observables p s i2L and p S 22L since such strategy 
might reveal the subtlety of determining the corresponding critical theory. Here, p s ; with i G {1,2} 
and L refer to the spin stiffness in the i direction and the spatial box size used in the simulations, 
respectively. Remarkably, while the Monte Carlo data we obtain for all the models considered here, 
including the herringbone- and ladder-dimer models on the square lattice, seem to be compatible 
with the recently proposed scenario of an enhanced correction to scaling, our data are in consistence 
with the established numerical values for the critical exponents in the 0(3) universality class as 
well. To explain the results presented in this study, a more detailed numerical study, and possibly a 
better theoretical understanding, for the critical theories of the models investigated here is required. 
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I. INTRODUCTION 



While being well-studied and understood thoroughly, 
the dimerized quantum Heisenberg models with spatial 
anisotropy have triggered theoretical interests again re- 
cently flH12|. For example, the 3-d spatially anisotropic 



quantum Heisenberg model with a ladder dimerization 
pattern is used to demonstrate a universal behavior, 
which is argued to be relevant for understanding the ex- 
perimental results of the material TICUCI3 Further, 
the 2-d dimerized spin-1/2 Heisenberg model with a spa- 
tially staggered anisotropy is of particularly interesting 
because this model seems to establish an unconventional 
phase transition [14| . Specifically, although it is believed 
that the phase transition induced by dimerization for this 
model should be governed by the 0(3) universality class 
theoretically |15l - [l9| . a recent large scale Monte Carlo 
calculation obtains v = 0.689(5) and f3/v = 0.545(5), 
which are in contradiction to the established 0(3) re- 
sults v = 0.7112(5) and (i/u = 0.519(1) in the literature 
[20j . Here v and f3 are the critical exponents correspond- 
ing to the correlation length and the magnetization, re- 
spectively. In order to clarify this issue further, several 
efforts have been devoted to study the phase transition 
of this model induced by dimerization. For instance, 
an unconventional finite-size scaling is proposed in [2lj . 
Further, in [22| it is argued that, due to a cubic term, 
there is a large correction to scaling for this phase tran- 
sition which results in the u nexp ected v = 0.689(5) and 
P/v = 0.545(5) obtained in (yf. Later, a Monte Carlo 
study indeed provides strong evidence to support this 



scenario of an enhanced correction to scaling [21[ . In ad 



dition to the staggered-dimer spin-1/2 Heisenberg model 
on the square lattice, in [22[ it concludes as well that 



a similar model on the honeycomb lattice, which is de- 
picted in the bottom panel of figure l 1 , as well as the 
herringbonc-dimer model on the square lattice (middle 
panel of figure 1) also belong to the same category of 
models receiving an enhanced correction. This general 
picture regarding the correction to scaling for 2-d dimer- 
ized quantum Heisenberg models is indeed supported by 
several related Monte Carlo studies H, I23T - I26T ] . 

While the Monte Carlo data of the staggered-dimer 
model on the square lattice provides convincing evidence 
for the proposal of an enhanced correction to scaling, 
the most noticeable observation regarding the critical 
behavior of this model is the good scaling property of 
the observable p S 2^L [2~l| . where again p S 2 and L are 
the spin stiffness in the 2-direction and the spatial box 
size employed in the simulations, respectively Inspired 
by this observation, one naturally would like to exam- 
ine whether for the staggcrcd-dimer model on the hon- 
eycomb lattice and the herringbone-dimer model on the 
square lattice, a similar good scaling behavior will be ob- 
served when considering the same observable p S 2 2£ for 
these two dimerized models. Although the phase tran- 
sition induced by dimerization of the staggered-dimer 
model on the honeycomb lattice has been studied be- 
fore, a detailed comparison between the scaling behavior 
of p s \2L and p S 2^L as well as the relevant investigation 
of the exponent /3/v are not available yet 2 . Further, to 
examine whether the enhanced correction to scaling, as 
suggested in (22j , does have impact on the determination 
of the exponents v and fi/v for the herringbone-dimer 
model is an interesting topic to explore as well 3 . Hence 
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1 We will call this model the staggered-dimer model as well unless 
confusion arises. 

2 An unpublished work of determining f}/v is available in [27f . 

3 An unpublished work of studying the critical theory for the 
herringbone-dimer model is available in [27I . 
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in this study, we have investigated the phase transitions 
of the herringbone- and ladder-dimer spin- 1/2 Heisen- 
berg models on the square lattice, as well as the quantum 
staggered-dimer model on the honeycomb lattice. In par- 
ticular, the largest lattice sizes reached here are as twice 
large as those of the relevant early studies in some cases. 
The results for the ladder-dimer model are included here 
for completeness and comparison purpose, since the en- 
hanced correction to scaling should be absent for this 
model. Remarkably, as we will demonstrate later, indeed 
p S 2^L of the staggered-dimer model on the honeycomb 
lattice shows a good scaling behavior. Consequently, we 
are able to obtain a value for v, in agreement with the 
established v = 0.7112(5) in the 0(3) universality class, 
by employing the leading scaling ansatz in our finite-size 
scaling analysis for p S 22L. Interestingly, while the Monte 
Carlo data we obtain for all the models studied here, in- 
cluding the herringbone- and ladder-dimer models on the 
square lattice, seem to be compatible with the recently 
proposed scenario of an enhanced correction to scaling for 
the phase transitions considered here, our data are also in 
consistence with the established results of v = 0.7112(5), 
fi/v = 0.519(1), and u> ~ 0.78 (a; is the confluent expo- 
nent) in the 0(3) universality class. To explain the re- 
sults presented in this study, a more detailed numerical 
study, and possibly a better theoretical understanding, 
for the critical theories of the models investigated here is 
required. 

This paper is organized as follows. First, after an in- 
troduction, the spatially anisotropic quantum Hciscnberg 
models and the relevant observables studied in this work 
are briefly described, after which we present our numer- 
ical results. In particular, the results obtained from the 
finite-size scaling analysis are discussed in detail. A final 
section then concludes our study. 

II. MICROSCOPIC MODEL AND 
CORRESPONDING OBSERVABLES 

The Hciscnberg models considered in this study are 
defined by the Hamilton operator 

H = ]T ■ S y + ^ J' S X > ■ Sy,, (1) 

{xy) {x'y') 

where J and J' are antiferromagnetic exchange couplings 
connecting nearest neighbor spins (xy) and (x'y'), respec- 
tively. Figure 1 illustrates the models which are described 
by Eq. ([I]) and are investigated in great detail here. To 
study the critical behavior of these models near the tran- 
sition driven by the anisotropy, in particular, to deter- 
mine the critical points as well as the critical exponent 
v, the spin stiffnesses in the 1- and 2-dircctions which are 
defined by 

Psi = ^{W?), (2) 

are measured in our simulations. Here /3 is the inverse 
temperature and L again refers to the spatial box size. 



Further (Wf) with i G {1,2} is the winding number 
squared in the i direction. In addition, the second Binder 
ratio Qi, which is defined by 



is also measured in our simulations as well. Here m z is 
the z component of the staggered magnetization rh s = 
J2x(~^) Xl+X2 S x - By carefully investigating the spa- 
tial volume and the J' /J dependence of p S iL as well as 
Q2, one can determine the critical points and the criti- 
cal exponent v with high precision. Finally the exponent 
P/v is determined by studying the scaling behavior of 
the observables (\m z s |) and ((m z s ) 2 ), which are measured 
in this study as well, at the corresponding critical points. 

III. DETERMINATION OF THE CRITICAL 
POINTS AND THE CRITICAL EXPONENT v 

To study the quantum phase transitions of our cen- 
tral interest, we have carried out lar ge s cale Monte Carlo 
simulations using a loop algorithm |28l432j |. Further, to 
calculate the relevant critical exponent v and to deter- 
mine the location of the critical points in the parameter 
space J' j J for the models described by figure 1, we have 
employed the technique of finite-size scaling for certain 
observables. For example, if a transition is second or- 
der, then near the transition the observable p s {lL for 
i € {1,2} and Q2 should be described well by the follow- 
ing finite-size scaling ansatz [33T - [37| 

L (t) = goitL 1 ^ ,L Z I (3,r) + L'^goJtL 1 ^ ,L Z / (3,r) 
= go {tLV»,L z /p,r) x 

(l + L-^g'oJtL^^^/fS,^) (4) 

where Ol stands for Q2 and p S iL with i £ {1,2}, L is 
the lattice size in the 1-direction, t = (j c — j)/j c with 
j = (J'/J), v is the critical exponent corresponding to 
the correlation length £, u> is the confluent correction ex- 
ponent, z is the dynamical critical exponent which is 1 
for the phase transitions considered here, and r is the 
ratio of the lattice size in the 1- and 2-direction. Fur- 
ther, go, go^ , and g' ^ appearing above are smooth func- 
tions of the variables tL 1 /" ' , L//3, and r. In practice one 
would carry out the analysis close to the critical point so 
that g' in Eq. (J4j) can be approximated by a constant. 
Specifically, the following ansatz 

L (t) = (1 + bL-^goitL 1 ^, L z /f3, r), (5) 

where b is some constant, is frequently used when apply- 
ing the finite-size scaling technique. While Eq. ([5]) is only 
valid for large box sizes and close to the critical point, to 
present the main results of this study we find it is suf- 
ficient to employ Eq. (J5j) for the data analysis. Notice 
that for square lattice or rectangular-shape lattice with 



3 



a fixed r, one will intuitively neglect the effect of r in 
Eq. ([5]). Hence, we will apply Eq. ([5]) with a constant r 
to the relevant observables for obtaining (J'/J) c and i/. 
Notice from Eq. ([5]), one concludes that the curves for 
Ol corresponding to different L, as functions of J'/ J, 
should intersect at the critical point (J'/J) c for large 
L. Without loss of generality, we have hxed J = 1 in 
our simulations and have varied J'. Additionally, the 
box size used in the simulations ranges from L = 24 to 
L = 136 (Strictly speaking, L = \/N for the staggered- J 
dimer model on the honeycomb lattice. Here N is the 
number of spins used in the simulations). To reach a lat- 
tice size as large as possible, we use /3J = 2L for each 
L in our simulations. Consequently, the temperature de- 
pendence in Eq. ((5]) drops out. First of all, let us focus on 
our results for the staggercd-dimer spin-1/2 Heiscnberg 
model on the honeycomb lattice. 

J' 

A. Results for the staggered-dimer spin-1/2 
Heisenberg model on the honeycomb lattice 



Figure 2 shows the Monte Carlo data of p s i2L, p S 2^L, 
and Qi with 24 < L < 96 as functions of J' /J for 
the staggered-dimer spin-1/2 Heiscnberg model on the 
honeycomb lattice. The figure clearly indicates that the 
phase transition is most likely second order since for all 
the observables p s i2L, p S 22L, and Q2, the curves of dif- 
ferent L tend to intersect near a particular point in the 
parameter space J'/ J. A surprising observation from fig- 
ure 2 is that, while p s \2L receives a sizable correction to 
its scaling (which has already been shown in [Ifjj]), the 
observable p S 2%L shows a good scaling behavior. Specif- 
ically, the correction to scaling for p S 2%L is negligible 
for L > 32. These findings are similar to the scenario 
regarding the correction to scaling for the same observ- 
ables, namely p a \2L and p S 22L of the staggered-dimer 



model on the square lattice |21|. Indeed, from p^L 
with L > 32, we are able to reach a value for v com- 
patible with the expected v = 0.7112(5) using the lead- 
ing finite-size scaling ansatz in Eq. ([5]) (letting b = in 
Eq. ([S])). For example, the v obtained from applying a 
second order Taylor expansion in tL l l v of Eq. ([5|) with 
b = to the observable p S 2%L with 40 < L < 96 is 
given by v = 0.7167(40), which is in nice agreement with 
its theoretical expectation in the literature. To reach a 
value for v consistent with v = 0.7112(5) using the lead- 
ing finite-size scaling ansatz and the observable p s i2L, 
one has to use data with fairly large L as indicated in 
[2^ | . Indeed a similar conclusion is reached here. Inter- 
estingly, with the observable p s \2L, while we either arrive 
at values of v statistically different from v = 0.7112(5) 
or cannot reach good results when b and w in Eq. ([5]) 
are included as fitting parameters, compatible results of 
v with v = 0.7112(5) can be obtained from the fits with 
the assumption that tu < 0.5 is used as a criterion for the 
fits. For instance, using p s \1L with 40 < L < 96, the 
values of v and u determined from the fits by cmploy- 




FIG. 1: The dimerized quantum Heisenberg models with spa- 
tial anisotropy considered in this study. 



ing the criterion of u> < 0.5 are given by v = 0.7054(45) 
and lu = 0.42(8), respectively. Notice u> ~ 0.42 is smaller 
than the expected uj ~ 0.78 in the 0(3) universality class. 
One might conclude that our results are consistent with 
the scenario outlined in [22[ that the correction to scal- 
ing for this model is enhanced due to a cubic irrelevant 
term. Finally, from Q2 with 32 < L < 96, a fit us- 
ing Taylor expansion to second order in tL 1 ^ as well as 
letting b = in Eq. © leads to v = 0.7102(56) and 
(J' /J) c = 1.73560(4), both of which agree quantitatively 
with the known results in the literature. Interestingly, 
the obtained coefficients for (tL 1 ^) 2 in the fits associ- 
ated with Q2 are very small. Hence, we can even reach a 
value for v in agreement with v = 0.7112(5) using a first 
order Taylor expansion in tL x l v of Eq. ([5]) to fit the data 
points of Q2 (the last 2 rows in tabic 1). The values of v 
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FIG. 2: Monte Carlo data of p s i2L, p s2 2L, and Q 2 with 24 < 
L < 96 for the staggered-dimer spin- 1/2 Heisenberg model on 
the honeycomb lattice. 

and (J'/J) c obtained from the fits mentioned above are 
listed in table 1. Notice that the uncertainties of (J'/J) c 
and v shown in table 1, as well as in tables 2, 3, 4, and 5 
in the following sections, are determined by a conserva- 
tive estimate based on the standard deviations obtained 
from the bootstrap resampling method employed for the 
fits. 



B. Results for the herringbone-dimer spin- 1/2 
Heisenberg model on the square lattice 

After having calculated (J'/J) c and v for the phase 
transition induced by dimerization of the staggered- 
dimer spin- 1/2 Heisenberg model on the honeycomb lat- 
tice, we turn to investigating the corresponding critical 
theory of the herringbone-dimer model on the square lat- 



observable 






L 




V {J'/J)c 


x 2 /dof 


p„i2L 


40 


< 


L 


< 


96 0.7054(45) 1.7355(2) 


1.2 


p a \2L 


48 


< 


L 


< 


96 0.7096(55) 1.7355(3) 


1.2 


Ps-2.1L 


32 


< 


L 


< 


96 0.7156(40) 1.73545(3) 


1.9 


ps2*2L 


40 


< 


L 


< 


96 0.7167(40) 1.73548(3) 


1.4 


Ps2lL 


48 


< 


L 


< 


96 0.7055(42)* 1.73550(3)* 


0.9 


ps2*2L 


56 


< 


L 


< 


96 0.7082(45)* 1.73551(3)* 


0.9 


Q2 


32 


< 


L 


< 


96 0.7102(56) 1.73560(4) 


1.4 


Q 2 


40 


< 


L 


< 


96 0.712(6) 1.73570(5) 


1.3 


Q2 


32 


< 


L 


< 


96 0.7107(56)° 1.73564(3)° 


1.4 


Q2 


40 


< 


L 


< 


96 0.712(6)° 1.73566(4)° 


1.4 



TABLE I: The numerical values of v and (J'/J) c calculated 
from p s \2L, p S 22L, and Q2 for the staggered-dimer model 
on the honeycomb lattice. All results are obtained by using 
a second order Taylor expansion in tL 1 ^ of Eq. ([5]) except 
those with a star (diamond) which are determined by a third 
order (first order) Taylor expansion. The confluent correction 
is included in the fit explicitly only for p s \2L and is assumed 
to satisfy the condition u> < 0.5. 



observable L v {J'/J)c xVDOF 

p s L 24 < L < 136 0.705(2) 2.49804(8) 1.4 

p s L 24 < L < 96 0.706(3) 2.4980(2) 1.4 

p s L 24 < L < 72 0.702(5) 2.4980(5) 1.1 

p s L 32 < L < 136 0.706(2) 2.49805(10) 1.4 

p s L 32 < L < 96 0.707(3) 2.4980(3) 1.4 

p s L 24 < L < 72 0.700(4)* 2.49813(10)* 1.1 

p s L 32 < L < 136 0.706(2)* 2.49806(3)* 1.5 

p s L 32 < L < 96 0.707(3)* 2.49806(7)* 1.4 

Q 2 24 < L < 136 0.714(4) 2.49800(7) 0.9 

Q 2 32 < L < 136 0.715(6) 2.4980(2) 1.1 

Q 2 24 < L < 136 0.710(4)* 2.49820(6)* 1.1 

Q 2 24 < L < 96 0.716(5)* 2.4983(1)* 1.1 

p s L 24 < L < 136 0.701(2) 2.49803(7) 1.2 

p s L 32 < L < 136 0.702(3) 2.4980(1) 1.2 

Q 2 24 < L < 136 0.707(4) 2.49800(15) 0.8 

Q2 32 < L < 136 0.7075(45) 2.4980(2) 0.9 

Q 2 32 < L < 136 0.7056(40)* 2.49810(6)* 0.9 

Q 2 40 < L < 136 0.7065(40)* 2.49810(7)* 0.8 



TABLE II: The numerical values of v and (J'/J) c calculated 
from p s L and Q2 for the herringbone-dimer model on the 
square lattice. While the results presented in the first twelve 
rows are obtained by using a second order Taylor expansion 
in tL 1 /" of Eq. ([5]), those listed in the last six rows are de- 
termined with a third order Taylor expansion. Further, all 
results are calculated with the ui and b in Eq. ([5} left as fit- 
ting parameters except those with a star which are determined 
through fits with a fixed uj = 0.78. 
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observable 






L 




V 


10 


x 2 /DOF 


p s L 


24 


< 


L 


< 


136 0.708(5) 


0.53(3) 


1.6 


p s L 


32 


< 


L 


< 


136 0.710(5) 


0.57(5) 


1.7 


p s L 


40 


< 


L 


< 


136 0.711(6) 


0.63(8) 


1.7 


PsL 


24 


< 


L 


< 


96 0.706(8) 


0.48(5) 


1.6 


p a L 


24 


< 


L 


< 


136 0.711(7)* 


N/A 


2.0 


p s L 


32 


< 


L 


< 


136 0.711(6)* 


N/A 


1.8 


Q2 


24 


< 


L 


< 


136 0.710(5) 


2.2(1) 


0.9 


Q2 


32 


< 


L 


< 


136 0.710(5) 


2.42(25) 


0.9 


Q2 


24 


< 


L 


< 


96 0.712(7) 


2.10(13) 


1.0 


Q2 


24 


< 


L 


< 


136 0.714(9)* 


N/A 


1.8 


Q2 


32 


< 


L 


< 


136 0.712(7)* 


N/A 


1.2 



TABLE III: The numerical values of v and (J'/J) c calculated 
from psL and Q2 for the herringbone-dimer model on the 
square lattice. All results are obtained by using a first order 
Taylor expansion in tL x l" of Eq. ((5} with cu and b left as fitting 
parameters except those with a star which are determined 
through fits with a fixed uj = 0.78. 



tice. Since for this model one has p s \ = p S 2, the rele- 
vant observables used in our finite-size scaling analysis 
arc p s L, which is the average of p s iL and p S 2L, and the 
second Binder ratio Q2 (figure 3). To calculate v, we first 
carry out several analysis by employing the second order 
Taylor expansion in tL 1 '" of Eq. ([5]), with the sublcading 
correction included explicitly, to fit our Monte Carlo data 
of p s L with variant range of L. Remarkably, a numerical 
value of v compatible with v = 0.7112(5) can be obtained 
if the smallest and largest box sizes used in the fits are 
larger than 24 and 96, respectively. The results of (J 7 / J) c 
and v calculated from these fits are listed as the first 5 
rows in table 2. Further, the values of uj determined 
from these fits ranges from 0.58 to 0.79 with an average 
of 0.66. Notice uj ~ 0.66 we obtain is slightly bellow the 
expected uj ~ 0.78 in the 0(3) universality class, hence is 
consistent with the scenario suggested in 22]. However, 
these results for uj should only be considered as effec- 
tive ones. Similarly, a fit using Q2 with 24 < L < 136 
as well as a second order Taylor expansion in tL 1 '" of 
Eq. ([5]), with the confluent correction left as fitting pa- 
rameters for the fit, leads to (J'/J) c = 2.49800(7) and 
v = 0.714(4). Notice the determined v = 0.714(4) is con- 
sistent with v = 0.7112(5). Further, the confluent expo- 
nent uj from the fit is given by uj = 2.0(2). Finally, while 
using other range of L for Q2 we can arrive at values of v 
agreeing with v = 0.7112(5), the uj calculated from these 
additional fits are poor determined. Notice that a sec- 
ond order Taylor expansion in tL x l v of Eq. (|5|) with the 
subleading correction included for the fit contains seven 
fitting parameters, which is at the border of reasonable 
amount of the unknown coefficients for a fit. Still, one 
would like to understand whether a consistent v with 
v = 0.7112(5) can be obtained from the fits with fewer 
fitting parameters. Interestingly, using (J'/J) c = 2.4980, 



the data points very close to the critical point, as well as a 
first order Taylor expansion in tL 1 ^ of Eq. ([5]) with the 
confluent correction included explicitly (which has five 
unknown coefficients), the values of v determined from 
these new fits for both p s L and Q2 are compatible with 
v = 0.7112(5) (table 3 and top panel of figure 4). There- 
fore we conclude that, our data points of p s L and Q2 for 
the herringbonc-dimcr model on the square lattice indeed 
can be described nicely with the expected v = 0.7112(5) 
in the 0(3) universality class. Notice that the values of uj 
calculated from the additional fits (first order Taylor ex- 
pansion of Eq. ([5])) related to p s L has an average of 0.55, 
hence again is in agreement with the scenario of a large 
correction to scaling for this phase transition induced by 
spatial anisotropy. 

Interestingly, while the results we obtain so far are in 
consistence with the scenario of a large correction to scal- 
ing as suggested in [22[ , the numerical values of v deter- 
mined from the fits with a fixed uj = 0.78 are also compat- 
ible with v = 0.7112(5) for both p s L and Q (table 2 and 
table 3). For instance, a fit using a fixed uj = 0.78 to the 
observable p s L with 32 < L < 96 leads to v = 0.707(3), 
which is in nice agreement with the expected result of 
v = 0.7112(5). Further, we are also able to arrive at val- 
ues of v agreeing with v = 0.7112(5) using the first order 
Taylor expansion in tL 1 /" of Eq. ([5]) with a fixed uj = 0.78 
for the fits (table 3). These additional fits contains only 
four unknown coefficients. Hence, both the strategies of 
fixing uj to be 0.78 or letting it be a fitting parameter 
lead to results of v consistent with v = 0.7112(5). Notice 
in our earlier calculations using a second order Taylor 
expansion of the full ansatz Eq. (JS|), although the mean 
average of uj, determined from the fits with both the b 
and uj in Eq. ([5]) left as fitting parameters, is smaller 
than the expected 0.78 in most of the cases, the uncer- 
tainties for uj from these fits are large. Hence, for a spread 
range (01,02) of uj (i.e. uj G (01,02)), consistent v with 
v = 0.7112(5) is obtained from the fits using the chosen uj 
in (ai, 02). Therefore, it might be premature to conclude 
that the confluent exponent uj for this phase transition is 
indeed smaller than 0.78 just from what we have obtained 
so far. In addition, since the values of uj obtained from 
the fits might be contaminated by higher order terms, 
a more sophisticated determination of uj should be per- 
formed. 



C. Results for the ladder-dimer spin- 1/2 
Heisenberg model on the square lattice 

The final dimer model considered in this study is the 2- 
d quantum Heisenberg model on the square lattice with a 
ladder spatial anisotropy, which has been studied exten- 
sively in the literature. Here for completeness, we have 
re-investigated the phase transition induced by dimer- 
ization of this model. Intuitively, due to their similarity, 
one might expect that the good scaling behavior of the 
observable p S 22L, found for the staggered- dimer model 
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observable 






L 






V (J'/J)c 


x 2 /dof 


ps\i.L 


24 


< 


L 


< 


136 


0.707(2) 1.90955(9) 


1.1 


p a \1L 


24 


< 


L 


< 


96 


0.705(3) 1.9096(2) 


1.2 


Pa\1L 


24 


< 


L 


< 


72 


0.696(5) 1.9095(6) 


1.1 


p a \2L 


32 


< 


L 


< 


136 


0.708(2) 1.90956(12) 


1.1 


Ps\1L 


32 


< 


L 


< 


96 


0.706(3) 1.9097(2) 


1.1 


p a \2L 


24 


< 


L 


< 


136 0.706(2)* 1.90960(3)* 


1.2 


Ps\1L 


32 


< 


L 


< 


72 


0.698(5)* 1.90956(16)* 


1.1 


p a \2L 


32 


< 


L 


< 


96 


0.707(3)* 1.90961(7)* 


1.2 


psi2L 


24 


< 


L 


< 


136 


0.704(3) 1.9095(1) 


1.1 


psi2L 


32 


< 


L 


< 


136 


0.705(3) 1.90956(13) 


1.1 



TABLE IV: The numerical values of v and (J'/J) c calculated 
from p s i2L for the ladder-dimer model on the square lattice. 
While the results presented in the first eight rows are obtained 
by using a second order Taylor expansion in tL 1 ^ of Eq. J5J, 
those listed in the last two rows are determined with a third 
order Taylor expansion. Further, all results are calculated 
with the to and b in Eq. ([5} left as fitting parameters except 
those with a star which are determined through fits with a 
fixed ui = 0.78. 



FIG. 3: Monte Carlo data of p 8 L and Q 2 with 24 < L < 136 
for the herringbone-dimer spin- 1/2 Heisenberg model on the 
square lattice. 



on both the square and honeycomb lattices, will emerge 
again for the ladder-dimer model. Interestingly, the ef- 
fects of the correction to scaling for p s \2L and p S 2%L 
of this model are about the same qualitatively (figure 
4). This indicates the fundamental difference regard- 
ing the correction to scaling between the staggered- and 
ladder-dimer models as suggested in [22| . Similar to the 
analysis performed for obtaining v and (J'/J) c for the 
herringbone-dimer model, a second order Taylor expan- 
sion in tL x l v of Eq. ([5]), with the subleading correction in- 
cluded explicitly, is employed to fit the Monte Carlo data 
of p s \2L with variant range of L. The obtained (J'/J) c 
and v are in table 4. Interestingly, table 4 implies that 
a result for v compatible with v = 0.7112(5) can be ob- 
tained as well if the smallest and largest box sizes used in 
the fits are larger than 24 and 96, respectively. Further, 
although with considerable large uncertainties, most the 
values of u> determined from these fits are smaller than 
0.78 and arc compatible in magnitude with those deter- 
mined from the herringbone-dimer model. In addition, 
using (J'/J) c = 1.9095, the data points very close to the 
critical point, as well as a first order Taylor expansion 
in tL 1 ^ of the full ansatz Eq. (|5|), the values of v deter- 
mined from these new fits for p s i2L are compatible with 
v = 0.7112(5) (table 5 and bottom panel of figure 4). 
Therefore we conclude that, our data points of p s \lL for 
the ladder-dimer model on the square lattice indeed can 
be described nicely with the expected v = 0.7112(5) in 
the 0(3) universality class. Finally, similar to the results 




2.495 2.496 




1.906 1.907 1.908 1.909 1.91 1.911 1.912 1.913 

J7J 



FIG. 4: Fits of p s L (24 < L < 136) of the herringbone- 
dimer model (top panel) and p a \2L (24 < L < 136) of the 
ladder-dimer model (bottom panel) to the first order Taylor 
expansion in tL 1 ^ of the full ansatz Eq. (J5]). While the circles 
are the numerical Monte Carlo data from the simulations, the 
solid curves are obtained by using the results from the fits. 



7 



2.6 - 




1.6 - 



1.904 1.906 1.908 1.91 1.912 1.914 




0.8 - 



1.906 1.908 1.91 1.912 1.914 

J7J 



FIG. 5: Monte Carlo data of p a \2L and p s2 2L with 24 < L < 
136 for the ladder-dimer spin-1/2 Heisenberg model on the 
square lattice. 



observable 






L 






x 2 /dof 


p 3 i2L 


24 


< 


L 


< 


136 0.707(6) 0.29(4) 


1.5 


p s \2L 


32 


< 


L 


< 


136 0.709(6) 0.26(6) 


1.5 


p 3 i2L 


24 


< 


L 


< 


96 0.703(8) 0.28(6) 


1.5 



TABLE V: The numerical values of v and (J'/J) c calculated 
from p a \2L for the ladder-dimer model on the square lattice. 
All results are obtained by using a first order Taylor expansion 
in tL x l" of the full ansatz Eq. 



for the herringbonc-dimer model, with a fixed u = 0.78, 
the values of v calculated from our finite-size scaling anal- 
ysis are compatible with v = 0.7112(5) as well (table 4). 



IV. DETERMINATION OF THE EXPONENT 

P/v 

After having calculated the critical exponent v from 
the relevant observables for the models described by fig- 
ure 1, we turn to the determination of the exponent (3 /v. 
To calculate P/v, the scaling behavior of the observables 
(|m*|) and ((mf) 2 ) are studied. Specifically, at critical 
points and for large L, the observable (\m z s \ k ) should scale 



as 

(\mt\ k ) = (a k + b k L-")L- k ^, (6) 

where a k , b k are some constants for k = 1 and each even 
positive integer k. Since precise knowledge of the criti- 
cal points is essential in determining the exponent P/v, 
we use the values of (J/J') c obtained in previous sec- 
tions when calculating the critical exponent v. Interest- 
ingly, as shown in tables 1, 2, and 4, small statistical 
deviation between some of the determined critical points 
of the same model is found. We attribute such small 
discrepancy to the presence of higher order sublcading 
corrections which are not taken into account in our anal- 
ysis, as well as the fact that the bootstrap resampling 
method used in calculating (J/J') c and v might occa- 
sionally leads to underestimated errors. While small de- 
viation is observed, the accuracy of (J/J') c presented in 
tables 1, 2, and 4, is sufficient for determining P/v by 
investigating the scaling behavior of (|mf |) and ((mf) 2 ), 
at the corresponding critical points. Hence, the values of 
( J I J') c for the herringbone- and ladder-dimer models on 
the square lattice, as well as the staggered-dimer model 
on the honeycomb lattice are taken to be 2.4980, 1.9095, 
and 1.7355, respectively. Further, we have carried out 
additional simulations at these critical points so that the 
largest lattice size we reach for both the herringbone- 
and staggered-dimer models is L = 184. First of all, let 
us focus on the results of P/v obtained from (|mf|). In- 
terestingly, with the expected leading scaling behavior, 
only from the ladder-dimer model we are able to reach a 
value of P/v which is in agreement with the established 
result P/v = 0.519(1) in the literature. For example, 
while a fit using the leading scaling expectation and the 
observable with L > 72 of the ladder-dimer model 

results in P/v = 0.517(2) (top panel of figure 6), the cor- 
responding numerical value of P/v determined from the 
same observable, with a similar range of L, is given by 
P/v = 0.527(3) (P/v = 0.531(3)) for the herringbonc- 
dimer model (staggcrcd-dimer model on the honeycomb 
lattice). Further, using L > 128 (L > 120), the value of 
P/u obtained from a fit without the confluent correction 
for the herringbone-dimer model (staggered-dimer model 
on the honeycomb lattice) is given by p/v = 0.522(5) 
(P/v = 0.526(3)). At this stage, one might conclude 
that the correction to scaling for the staggered- and 
herringbone-dimer models are indeed enhanced as pro- 
posed in (22j . In particular, the value of 10 determined 
from the related fits should be smaller than the ex- 
pected 0(3) result tu ~ 0.78. Surprisingly, for both 
the herringbonc-dimer model on the square lattice and 
the staggered-dimer model on the honeycomb lattice, us- 
ing the data of (|mf|) with 16 < L < 184, a fit in- 
cluding the subleading correction and a fixed u> = 0.78 
leads to values of p/u compatible with P/v = 0.519(1) 
(middle and bottom panels of figure 6). Further, if uj 
is left as a fitting parameter, although consistent P/v 
with P/v = 0.519(1) are obtained from the fits associ- 
ated with the herringbonc-dimer model, the uncertainties 
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for p/v and uj arc increased significantly. For example, 
the P/v and uj determined from a fit using (|mf|) with 
16 < L < 184 of the herringbonc-dimer model are given 
by p/v = 0.521(9) and 0.86(50), respectively. Finally, 
a fit to the observable (\m z s |) of the ladder-dimer model 
with a fixed ui = 0.78 leads to fj/v = 0.516(3), which 
is consistent with P/v = 0.519(1) as well. Interestingly, 
a value of P/v slightly below P/v = 0.519(1) is reached 
when the leading scaling prediction is employed to fit 
all available data of (Im^l) of the ladder-dimer model. 
This in turn implies that the coefficient b\ in Eq. ([5]) 
for the ladder-dimer model is small in magnitude. In- 
deed, the magnitude of b\ obtained from applying the 
full ansatz Eq. j6]) with a fixed ui = 0.78 to (|mf|) of the 
ladder-dimer model is of order 10~ 2 (The uncertainty for 
&i is comparable to b\ in magnitude as well). Finally, we 
have also carried out an additional analysis with different 
fixed values of uj in the fits. The obtained P/v for these 
additional fits are shown in table 5. From table 5 one 
concludes that the values of U) that would lead to consis- 
tent P/v with P/v = 0.519(1) for both the herringbonc- 
and staggcrcd-dimcr models ranges from 0.7 to 0.9, which 
matches reasonable well with the expected value 0.78. 
All the results we have reached so far imply that, our 
data points of (|m^|) for all the three dimerized models 
depicted in figure 1 are compatible with the established 
value of uj ~ 0.78 in the 0(3) universality class. 

After having demonstrated that our Monte Carlo data 
of (|mf|), for all the three dimerized models investi- 
gated here, are compatible with the established results 
of P/v = 0.519(1) and ui = 0.78 in the 0(3) universality 
class, a similar scenario is reached when considering the 
observable ((mf) 2 ). For example, using the leading scal- 
ing prediction, only from the ladder-dimer model we can 
reach a value for p/v compatible with P/v = 0.519(1) 
(top panel of figure 7). Further, with a fixed uj = 0.78, 
the fits for the herringbone-dimer model on the square 
lattice and the staggered-dimer model on the honeycomb 
lattice result in P/v = 0.5202(15) and P/v = 0.518(2), 
respectively (middle and bottom panels of figure 7) . The 
results of p/v = 0.5202(15) and p/v = 0.518(2) we just 
obtain arc in quantitative agreement with the expected 
P/v = 0.519(1). Further, an analysis for ((mf) 2 ), with 
variant fixed values of uj in the fits, leads to a similar 
conclusion like that of (|mf |). Specifically, from ((ml) 2 }, 
the values of uj that would lead to consistent p/v with 
P/v = 0.519(1) for both the herringbone- and staggered- 
dimer models ranges from 0.7 to 0.9 as well (table 6). 

Tables 4, 5, and 6 summarizes our calculations on 
determining p/v for the phase transitions induced by 
dimcrization for all the three models shown in figure 1. 



V. DISCUSSION AND CONCLUSION 

In this study, we investigate the phase transitions in- 
duced by dimcrization for the herringbone- and ladder- 
dimer spin-1/2 Heiscnbcrg model on the square lattice, 
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FIG. 6: Determination of p/v for the ladder- and 
herringbone-dimer models (top and middle panels) on the 
square lattice as well as for the staggered-dimer model on the 
honeycomb lattice (bottom panel). While P/v = 0.517(2) for 
the ladder-dimer model is obtained by fitting the 7 largest L 
data points of (|mf |) to their expected leading scaling behav- 
ior, fits with a fixed ui — 0.78 to |) of the herringbone- 
dimer model on the square lattice and the staggered-dimer 
model on the honeycomb lattice result in j3/v = 0.520(2) and 
P/v = 0.518(2), respectively. 



as well as the staggered-dimer model on the honeycomb 
lattice. In particular, we determine the values of the 
exponents v and P/v with high accuracy by employing 
the finite-size scaling analysis to the relevant observables. 
Similar to the scenario found for the staggered-dimer 
model on the square lattice, while the observable p s ±2L 
of the staggered-dimer model on the honeycomb lattice 
receives a sizable correction to its scaling, the observable 
p S 2^L shows a good scaling behavior. As a result, us- 
ing the data points of p a 2^L with moderate lattice sizes 
as well as the corresponding leading finite-size scaling 
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0.021 
0.018 
A 0.015 
oi 0.012 

E 

y 0.009 
0.006 
0.003 



.E. 
V 



0.018 
0.015 
A 0.012 

M 

M ^<» 0.009 
V 0.006 



(Vv = 0.515(2) 
X 2 /DOF = 0.45 



p/v = 0.5202(15) 
X 2 /DOF = 1 .0 
co = 0.78 (fixed) 



30 60 90 120 150 180 



p/v = 0.518(2) 
X 2 /DOF = 1 .4 
co = 0.78 (fixed) 



FIG. 7: Determination of Pjv for the ladder- and 
herringbone-dimer models (top and middle panels) on the 
square lattice as well as for the staggered-dimer model on the 
honeycomb lattice (bottom panel). While pfv = 0.515(2) for 
the ladder-dimer model is obtained by fitting the 5 largest L 
data points of ((mf ) 2 ) to their expected leading scaling behav- 
ior, fits with a fixed ui = 0.78 to ((m z a ) 2 ) of the herringbone- 
dimer model on the square lattice and the staggered-dimer 
model on the honeycomb lattice result in Pfv = 0.5202(15) 
and Pjv = 0.518(2), respectively. 



model 


L 


P/v 


X 2 /DOF 


ladder 


72 < L < 136 


0.517(2)* 


0.75 


ladder 


16 < L < 136 


0.516(3) 


1.1 


herringbone 


72 < L < 136 


0.527(3)* 


1.5 


herringbone 


128 < L < 184 


0.522(5)* 


0.7 


herringbone 


16 < L < 184 


0.520(2) 


1.0 


staggered 


72 < L < 136 


0.531(3)* 


1.2 


staggered 


120 < L < 184 


0.526(3)* 


0.55 


staggered 


16 < L < 184 


0.518(2) 


1.35 


ladder 


88 < L < 136 


0.515(2)* 


0.45 


ladder 


16 < L < 136 


0.515(2) 


1.3 


herringbone 


128 < L < 184 


0.525(4)* 


1.25 


herringbone 


16 < L < 184 


0.5202(15) 


1.0 


staggered 


128 < L < 184 


0.527(4)* 


1.1 


staggered 


16 < L < 184 


0.518(2) 


1.4 



TABLE VI: The numerical values of P/v calculated from 
(|mf|) (the first 8 rows) and ((mf) 2 ) (the last 6 rows) for 
the dimerized models considered in this study. All results 
are obtained with a fixed uj = 0.78 except those with a star 
which are determined by using the expected leading scaling 
prediction. 



model 


u (fixed) 


p/u 


X 2 /DOF 


herringbone 


1.0 


0.523(2) 


1.0 


herringbone 


0.9 


0.522(2) 


1.0 


herringbone 


0.7 


0.5181(23) 


1.06 


herringbone 


0.65 


0.5168(25) 


1.05 


herringbone 


0.6 


0.5152(27) 


1.04 


herringbone 


0.45 


0.5078(36) 


1.05 


staggered 


1.0 


0.523(2) 


1.55 


staggered 


0.9 


0.521(2) 


1.46 


staggered 


0.7 


0.5150(25) 


1.3 


staggered 


0.65 


0.5127(25) 


1.23 


staggered 


0.6 


0.5101(25) 


1.2 


staggered 


0.45 


0.4974(35) 


1.15 



TABLE VII: The numerical values of Pjv calculated from 
(|mf|) (16 < L < 184) with variant fixed ui for the fits. 



ansatz (letting b = in Eq. ([5])), we are able to obtain 
a value for v consistent with the expected v = 0.7112(5) 
in the 0(3) universality class. To understand this obser- 
vation for the staggered-dimer model on both the square 
and honeycomb lattices from field theory aspect, is an in- 
teresting and important topic to explore. In particular, 
whether the cubic term introduced in (22l | is responsible 
for this unexpected result should be investigated. Fur- 
ther, while it is argued in [22| that the herringbone-dimer 
model belongs to the category of models receiving a large 
correction, and our investigation supports this scenario, 
namely the values of ui obtained from the fits related to 



p s L are all smaller than the expected uj = 0.78, the same 
observable p s L as well as Q2 are also compatible with 
the established result of u> ~ 0.78 in the 0(3) universal- 
ity class. Specifically, with a fixed oj = 0.78, we are able 
to arrive at consistent v with v = 0.7112(5) from both 
p s L and Q2 of the herringbone-dimer model. In order 
to clarify whether our data is really compatible with the 
scenario of an enhanced correction to scaling suggested in 
[22^ , it will be desirable to carry out a more detailed in- 
vestigation to determine uj with high precision. In partic- 
ular, the consistence of the v, obtained from the fits using 
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model u(fixed) P/v xVDOF 



herrin 


gbone 


1.0 


0.5232(12) 


1.03 


herrin 


gbone 


0.9 


0.5221(13) 


1.0 


herrin 


gbone 


0.7 


0.5185(15) 


1.0 


herrin 


gbone 


0.65 


0.5172(16) 


1.0 


herrin 


gbone 


0.6 


0.5156(18) 


1.05 


herrin 


gbone 


0.45 


0.5082(25) 


1.1 


stagg 


ered 


1.0 


0.5235(15) 


1.7 


stagg 


ered 


0.9 


0.5213(15) 


1.55 


stagg 


ered 


0.7 


0.515(2) 


1.4 


stagg 


ered 


0.65 


0.512(2) 


1.35 


stagg 


ered 


0.6 


0.509(2) 


1.35 


stagg 


ered 


0.45 


0.4920(36) 


1.5 



TABLE VIII: The numerical values of ft/v calculated from 
{{m z a ) 2 ) (16 < L < 184) with variant fixed uj for the fits. 



els with spatial anisotropy considered here are compat- 
ible with the established results of f3/v = 0.519(1) and 
uj ~ 0.78 in the 0(3) universality class as well. Finally, 
the consistence of the (3/v, determined from the fits using 
a fixed uj G {0.7,0.9}, with fi/v = 0.519(1) implies that 
the observed enhanced correction to both the staggered- 
and herringbone-dimer models is because of the nonuni- 
versal coefficients bk in Eq. ©. Indeed, in tables 6, 7, and 
8, the values of b\ and 62 determined from the fits associ- 
ated with the herringbone- and staggered- dimer models 
are at least several times larger in magnitude than those 
of the ladder-dimer model. Whether there is a subtlety 
behind this observation or it is just a coincidence should 
be investigated since the determination of uj might be 
contaminated by higher order terms. In conclusion, to 
explain what have been obtained in this study, a more 
detailed numerical study, and possibly a better theoreti- 
cal understanding, for the critical theories of the models 
investigated here is required. 



a fixed uj = 0.78, with v = 0.7112(5) as shown in tables 2, 
3, 4 and 5, is unlikely a coincidence considering the fact 
that the conclusion is valid for both the observables spin 
stiffness and second Binder ratio of both the herringbonc- 
and ladder-dimer models. In addition, the finite-size scal- 
ing analysis performed for the determination of fi/v sug- 
gests that, our data points for all the 2-d dimerized mod- 
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